Quantitative rescattering theory for laser-induced high-energy plateau photoelectron 

spectra 
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A comprehensive quantitative rescattering (QRS) theory for describing the production of high- 
energy photoelectrons generated by intense laser pulses is presented. According to the QRS, the 
momentum distributions of these electrons can be expressed as the product of a returning electron 
wave packet with the elastic differential cross sections (DCS) between free electrons with the target 
ion. We show that the returning electron wave packets are determined mostly by the lasers only, 
and can be obtained from the strong field approximation. The validity of the QRS model is carefully 
examined by checking against accurate results from the solution of the time-dependent Schrodinger 
equation for atomic targets within the single active electron approximation. We further show that 
experimental photoelectron spectra for a wide range of laser intensity and wavelength can be ex- 
plained by the QRS theory, and that the DCS between electrons and target ions can be extracted 
from experimental photoelectron spectra. By generalizing the QRS theory to molecular targets, 
we discuss how few-cycle infrared lasers offer a promising tool for dynamic chemical imaging with 
temporal resolution of a few femtoseconds. 

PACS numbers: 32.80.Rm, 32.80.Fb, 42.50.Hz 



I. INTRODUCTION 

Much of our knowledge of the nonlinear interaction of 
intense laser radiation with atoms and molecules comes 
from the study of above-threshold ionization (ATI) which 
is characterized by a sequence of peaks in the electron 
spectrum, spaced by the photon energy. Since its first 
observation the subject has been "reinvestigated" 
many times. In 1987, as sub-picosecond laser pulses be- 
came available, it was shown 2] that ATI peaks suffer 
significant energy shifts and broadening, and each peak 
breaks up into substructures due to resonance enhance- 
ments produced by ponderomotive shifts of states. These 
substructures are called Freeman resonances. The nature 
of these Freeman resonances have been carefully investi- 
gated, for example for Ar in Wiehle et al. [3]. In re- 
cent years, with the introduction of COLTRIM detectors 
where electrons are measured over almost the whole 4n 
angular region, the two-dimensional (2D) electron mo- 
mentum spectra or the longitudinal electron momentum 
spectra of the photoelectrons or the target ions have been 
reported 0, H[ • These measurements reveal considerable 
structure not only in the electron energy distributions, 
but also in the angular distributions. All of these stud- 
ies focus on low energy electrons which are generated 
either by multiphoton ionization mechanism or by the 
tunneling ionization mechanism. According to the "con- 
ventional" wisdom, depending on the Keldysh parameter, 
7 = W Ip/2U P where I p is the ionization energy of the tar- 
get and Up the ponderomotive energy, if 7 is larger than 
one, the ATI electrons are generated by multiphoton pro- 



cesses, while if 7 is small, tunneling ionization is respon- 
sible for producing the low energy electrons. However, 
such a distinction is by no means clear-cut. In Ref. [f|, 
the 2D electron momentum spectra display pronounced 
fan-like structures even for laser intensities well into the 
tunneling ionization regime. Theoretical studies of elec- 
tron momentum spectra 0, 0] obtained from solving the 
time-dependent Schrodinger equation (TDSE) show that 
even in the tunneling region, photoelectron spectra show 
features that can be identified with the absorption of in- 
tegral number of photons, as revealed by the angular dis- 
tributions of the electron at fixed energies. These theo- 
retical multiphoton ionization features, when convoluted 
with the effect of the spatial distribution of intensities in 
a focused laser beam, can well reproduce the observed 
experimental 2D electron momentum spectra 0. 

While low-energy electrons, with energy less than 
about 2U P , account for the majority of the electrons gen- 
erated by an intense laser, already since 1993 photoelec- 
trons extending to 1QU P or more have also been observed. 
These electrons, unlike the low energy electrons gener- 
ated by the nonlinear processes, do not change with the 
electron energy rapidly, until a new cutoff near about 
lOUp is reached. They are known as high-energy plateau 
photoelectrons. Experiments showed that these electrons 
exhibitpronounced sidelobes not seen in low-energy elec- 
trons SEE 111 El- Th ese high-energy ATI (HATI) elec- 



trons have been interpreted as due to the rescattering 
process [3]. According to this model, electrons that are 
freed from the target atom at some well-defined ioniza- 
tion time may be driven back to revisit its parent ion. If 
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these returning electrons are backscattered by the target 
ion, they can be further accelerated by the laser field and 
emerge as high-energy electrons, reaching up to about 
lOUp. However, the plateau electron spectra, with ener- 
gies from AU P to about 10J7 P , are not always similar for 
different targets. For target like xenon, the plateau is 
flat, but for other like krypton, the plateau drops steeply 
as the electron energy increases. These features actu- 
ally change with peak laser intensities. Further studies 
of these HATI spectra around 1997 on inert gases discov- 
ered that resonantlike enhancements occur in the elec- 
tron spectra for particular laser intensities. Depending 
on the inert gas used, separate series of peaks have been 
observed. These observations have generated a flurry of 
theoretical interest Models based on analyzing re- 

sults from solving TDSE [H HH and from Floquet the- 



ory [17|, [18( have been proposed. Others are based on the 



channel-closing theory [19J, [20j] . More recent experiments 
confirmed that these enhancement disappears when laser 
pulse duration is reduced [1(3] • Since the HATI yields are 
four to five orders of magnitude smaller, despite of these 
experimental investigations, there are few systematic the- 
oretical calculations in the literature. In recent years, 
HATI electrons have drawn attention again since when 
they are generated by few-cycle laser pulses, their counts 
on the left and the right detectors along the laser polar- 
ization axis are different. Such asymmetry can be used 
to determine the absolute value of the carrier-envelope- 
phase (CEP) of the few-cycle pulses (2l| . 

Recently, we investigated the 2D high-energy photo- 
electron momentum spectra for atomic targets within 
the single active electron approximation based on the 
well-known rescattering model [H, [23| ■ We proposed a 



quantitative rescattering (QRS) theory 24, 25| where the 
HATI electrons are modeled as due to the backscattering 
of the returning electrons by the target ion. According 
to the QRS, high-energy photoelectron momentum dis- 
tributions D(k, 9) is shown to be expressed simply as 



D(k, 6) = W{k r )a{k ri 9 r ) 



(1) 



where a(k r ,9 r ) is the elastic differential cross sections 
(DCS) between free electrons, with momentum k r , with 
the target ion. Here 9 r is the scattering angle with re- 
spect to the the direction of the returning electrons along 
the laser polarization axis. In this equation, W(k r ) is in- 
terpreted as the momentum distribution of the returning 
electrons, to be called returning wave packet (RWP) in 
this paper. The validity of this QRS model has been 
tested using D(k,9) calculated from solving the TDSE, 
and a(k r ,9 r ) from the standard quantum mechanical 
scattering theory. Since there is a one-to-one relation be- 
tween (k, 9) and (k r , 9 r ), there are a number of important 
results from the QRS, as reported in our recent papers. It 
was shown in Morishita et al. 24 1 and Chen et al. [25( that 
one can extract elastic scattering cross sections a(k r ,9 r ) 
between free electrons and atomic ions from the HATI 



electron momentum spectra. The predictions have been 
confirmed in three recent experiments [1^, 27, [28[. In 
Chen et al. [29( ] . it was further shown that the momen- 
tum distribution W(k r ) of the RWP can be extracted 
from the second-order strong field approximation (SFA2), 
and that W(k r ) depends very little on the target (i.e., up 
to an overall normalization which is related to the total 
ionization probability). Thus one can use SFA2 to ob- 
tain W(k r ). By multiplying it with a(k r ,9 r ), we can 
use Eq. |T]) to obtain accurate high-energy photoelectron 
momentum distribution D(k,9). Since W(k r ) depends 
mostly on the lasers only, the target dependence of the 
HATI spectra can thus be explained based on the behav- 
ior of the elastic scattering cross sections a{k r ,9 r ). Based 
on the QRS model, the energy dependence of plateau ATI 
electrons seen for different targets are easily understood 
[25l ]. For a given target but different lasers, the HATI 
momentum spectra are determined by the RWP. Apply- 
ing the QRS model to few-cycle pulses, where the RWP 
varies with the change of the CEP, we have shown that 
the QRS theory can be easily used to retrieve the ab- 
solute value of the CEP [28|, Since nonsequential 
double ionization of atoms and molecules are understood 
based on the rescattering mechanism, in Micheau et al. 
[3lT | we showed that using the wave packet W{k r ) ob- 
tained from the HATI spectra, we can use QRS to obtain 
nonsequential double ionization yields. 

In this paper, we provide the full details of the QRS 
theory on atomic targets within the single active elec- 
tron approximation and establish its validity. Clearly our 
goal is not to limit ourselves to atomic targets only. We 
would like to study HATI spectra from molecular targets 
as well, in particular, from transient molecules. Recall 
that HATI electrons result from backscattering of the re- 
turning electrons by the target ion, i.e., electrons that 
undergo hard collisions with the target. Thus one should 
be able to retrieve the structure information of the target 
from the HATI spectra. In fact, we have shown that this 
is indeed possible for atomic targets already (H)]. Since 
laser pulses of duration of a few femtoseconds are read- 
ily available, one can perform pump-probe measurements 
where the pump pulse initiates a chemical reaction, such 
that the atomic coordinates of the molecule would evolve 
in time. Using a probe laser to take HATI spectra at 
different time delays, one will then have the opportu- 
nity to extract the structure of the transient molecule 
as a function of time from the measured HATI spectra. 
Thus short laser pulses may serve as a powerful tool for 
dynamic chemical imaging of transient molecules, with 
temporal resolution of a few femtoseconds. 

The rest of this paper is arranged as follows. In the 
next section, we discuss how to calculate the HATI elec- 
tron momentum spectra by solving the TDSE and us- 
ing the SFA2. We also explain how the elastic differ- 
ential cross sections are calculated. We then establish 
the QRS model. The validity of the QRS model is care- 
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fully examined in Section III, by testing against results 
obtained from solving the TDSE. In Section IV we illus- 
trate the application of the QRS model to experimental 
HATI spectra. To compare with experimental electron 
energy spectra, we include the laser focus volume effect. 
We finish the paper with a summary and outlook. We 
point out that a similar QRS model has been developed 



Born time zone Returning time zone 



for high-order harmonic generation (HHG) [24 . l33l . 

Atomic units are used in this paper unless otherwise 
noted. We also mention that in all the calculations the 
CEP is set to zero in this paper. 



II. THEORETICAL BACKGROUND 

The theory part is separated into four sections. Since 
the concept of rescattering can be understood in classical 
mechanics, we first consider the classical rescattering the- 
ory for an electron in a one-dimensional (ID) monochro- 
matic laser field. Then we discuss the calculations of ATI 
spectra by solving the TDSE and using the SFA2. For 
completeness we also include how the elastic scattering 
cross sections are computed. 

A. Classical one-dimensional rescattering theory 

A classical ID rescattering theory has been discussed 
by Paulus et al. 13]. Suppose that an electron in the ID 
atom is first released at some time to into a monochro- 
matic laser field E(t) = zEo cos cut, the Newton's equa- 
tion of motion for this system is given by 

z(t) = —Eq costot, (2) 

consequently the position of the electron at time t reads 

/ \ / \ -^0 r l 

z(t) = z(t ) H — [cos cot — cos (jjto\ 

LO A 



— - sin win + i(io) 

10 



(t-t ). 



(3) 



If the electron is initially at the origin with zero initial 
velocity, the time t r at which it returns to the origin 
satisfies 



costot r — coswtp + u> sin Lot Q (t r — t ) = 0. 



(4) 



The electron will never return to the origin if it is ionized 
before the laser field reaches its peak value while it can 
return to the origin more than once when it is born after 
the peak. It has been shown that higher order returns 
make very small contribution to the yield [29]. We con- 
sider only the first return here. The momentum k r of the 
electron when it first returns to the origin at time t r is 



E 

k r = z(t r ) = (sin ujt r + sinwto) ■ 

LO 
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FIG. 1: (Color online) Classical model of a ID electron in a 
monochromatic laser field, (a) Electric field and vector po- 
tential for a laser at the intensity of 1.0 x 10 14 W/cm 2 and 
wavelength of 800 nm. The born time and returning time 
zones are marked, (b) The electron velocity k r , vector po- 
tential A, and their ratio at the time of return t=t r , within 
the returning time zone, (c) Photoelectron energy after it has 
been backscattered by an angle 8 r = 180° against the return- 
ing time (bottom horizontal axis). Also shown is the relative 
ionization rate for electrons released with respect to the born 
time (top horizontal axis). 



In Fig. 1(a) the electric field and the vector potential 
are plotted vs Lot for a typical set of laser parameters, 
say I = 1.0 x 10 14 W/cm 2 and A = 800 nm. The max- 
imum value of the vector potential is Ao = 0.94. In the 
rescattering model, consider electrons that return within 
231° < iot r < 309°. In this region, A(t r ) > 0.78 A) and 
by solving Eq. (3), the corresponding born time is within 
4° < Loto < 25°. Both the born time zone and returning 
time zone are marked in Fig. 1(a). The electron born at 
time Loto = 13° returns at time tot r — 270°. It should be 
noted that electrons born before 13° return at the time 
after 270° and they follow a long trajectory, while those 
born after 13° return before 270° and they follow a short 
trajectory. 

The vector potential, the velocity of the electron at the 
return time, and the ratio of the returning velocity to the 
vector potential are plotted in Fig. 1(b), in the returning 
time zone. It can be seen that k r (t)/A(t) decrease from 
1.76 at 220° to 1.07 at 320°. 

We next assume that the electron undergoes collision 
with the target ion and elastically scattered by an angle 
8 r with respect to its incident direction. For t > t r , the 
components of its velocity, along the polarization axis 
and perpendicular to it, are given by 

z(t) = [smtot — smtoti + cosf3 r (smiot r sinwi )] , 

to 



(5) 



y(t) — (smiot r — sinwto) • (6) 

to 

From the above equation, the photoelectron energy Ek 
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measured by the detector outside of the field can be 
obtained by subtracting the ponderomotive potential, 
Up = Eq/(4uj 2 ), from the time averaged kinetic energy 

Ek = 2t/ p [sin 2 ujto + 2sbxu>t r (l — cos 9 r )(smu>ti — sinwio)]-( 

It is easy to find from ([7]) that, if the electron is born 
at uito = 14°, it returns to the origin at cot r = 265° 
when the vector potential almost reaches the maximum, 
then the photoelectron will have the maximal energy 
i?™ ax = 10.007L/p provided that the returning electron 
experiences a backwardscattering of 8 r — 180°. The elec- 
tron's kinetic energy Ek measured by the detector along 
the polarization axis is shown in in Fig. 1(c), together 
with the returning time versus with the relative ioniza- 
tion rate of electrons released in the born time zone. The 
lower horizontal axis in Fig. 1(c) indicates the returning 
time while the corresponding born time is given on the 
top. In the QRS model, we investigate backscattered 
electrons with energies greater than 4U p . From Fig. 1(c), 
it can be seen that electrons return at two different times 
could have the same kinetic energy except for the time 
around 265° at which Ek has the maximal value. Ac- 
tually, electrons which return after 265° are born before 
14°. These electrons have higher ionization rate than 
those born after 14°. For example, the ionization rate 
for electrons which are born at 22° and return around 
240° is about 3 times smaller than electrons born at 2° 
and return at 320°. 

Let us look back at Fig. 1(b) which shows that the ratio 
k r (t)/A(t) becomes flatter after a sharp decrease until 
240° and the mean value of k r (t)/A(t) in the returning 
time range from 240° to 320° is about 1.25. Although the 
ratio k r (t) /A(t) for returning time less than 240° deviates 
more from this mean value, the electron yield from this 
part is very small, as seen in Fig. 1(c). 



B. Method of solving the time-dependent 
Schrodinger equation 

The method for solving the time-dependent 
Schrodinger equation has been described in our 
previous works @, [|| where we studied the low-energy 
electron momentum spectra. Much more effort is needed 
to obtain accurate momentum spectra for high energy 
electrons. Here we describe the essential steps of the 
calculations. 

We treat the target atom in the single active electron 
model. The Hamiltonian for such an atom in the presence 
of a linearly polarized laser can be written as 



H = H + Hi(t) = -iv 2 + V(r) + H t {t). (8) 

The atomic model potential V(r) is parameterized in the 
form 



V{r) 



1 + a\e: 



+ a 3 re 



+ a 5 e 



(9) 



The parameters in Eq. ([9]) are obtained by fitting the 
calculated binding energies from this potential to the ex- 
perimental binding energies of the ground state and the 
jirst few excited states of the target atom. The param- 
eters for the targets used in this paper can be found in 
|35l ]. For Kr and Xe, we use the potential given by Gar- 
vey et al. 36]. The model potential for a neutral atom 
can also be expressed as 



V(r) = V s (r) - 1/r, 



(10) 



where V s (r) is a short-range potential. The atom-field 
interaction Hi(t), in length gauge, is given by 



ffi = r-E(t). 



(11) 



For a linearly polarized laser pulse (along the z axis) with 
carrier frequency u> and the CEP, ip, the field is taken to 
have the form 

E(t) = zE cos 2 ^— cos(ut + ip) (12) 

for the time interval (— r/2, r/2) and zero elsewhere. The 
pulse duration, defined as the full width at half maximum 
(FWHM) of the intensity, is given by T = r/2. 75. 

The time evolution of the electronic wavefunction 
*(r,t), which satisfies the TDSE, 

. d 



.^*(r,t) 



H^{r, t) 



(13) 



is solved by expanding in terms of eigenfunctions, 
R n i(r)Yi m (r), of H , within the box of r E [0, r max ] 



*(r,t) =J2c n t(t)R nl (r)Y lm {r) 



(14) 



ni 



where the radial functions R n i(r) are expanded by the 
DVR (discrete variable representation) [37, 38|, 3^] basis 
set associated with Legendre polynomials, while the c n i 
are calculated using the split-operator method [4(| 



c nl (t + At) ~ J2{e- lHoAt/2 e~ lH 



i(t+At/2)At 



xe 



iH At/2 



At) (15) 



where the matrix elements are evaluated efficiently by us- 
ing the DVR quadrature. For short pulses, say, T = 8 fs, 
and for electron energy as high as 12U P , converged re- 
sults can be obtained by setting r max = 1200. Note that 
in (|14|) . only m = is taken into account since for lin- 
early polarized laser pulses, contribution to the ioniza- 
tion probability from m — ±1 is relatively much smaller 
in comparison to the m = component. 

The photoelectron yield is computed at the end of the 
laser pulse by projecting the total final wave function 
onto eigenstates of a continuum electron with momentum 
k, 



D(k,0) 



Q3p 



|($j^(< = r/2))| 2 



(16) 



5 



where the continuum state satisfies the following 
equation 



V(r) 



4>, 



(17) 



C. Strong field approximation for calculating ATI 
electron spectra 

While direct solution of the Schrodinger equation in a 
time-dependent laser field has been widely used, the sim- 
pler strong field approximation is of interest for analyzing 
features of intense laser-atom interactions. 

By treating electron-laser interaction as the strong field 
and electron-target ion interaction as a perturbation, the 
amplitude for generating a photoelectron with momen- 
tum k is given by 



/(k)=/ 1 (k) + / 2 (k) 



(18) 



where we included the first two terms of the perturbation 
series only. 

In Eq. (|18p , the first order term /i (k) , which is tradi- 
tional called the strong field approximation (SFA), will 
be called SFA1, to be distinguished from the second term 
/ 2 (k), which will be called SFA2. The SFAl is given by 

/OO 
<ft(Xk(*)|ffi(*)l*o(*)) (19) 
-oo 

where is the ground state wavefunction. The Volkov 
state Xk is given by 



1 



(2tt) 3 / 2 
where the action S is 



i[k+A(t)].r -iS(k,t) 



(20) 



S(k,t) = ±J dt'[k + A(t')f. 
The SFA2 term in fTB|) is expressed as 

/oo poo p 

dt dt' dp( Xk (t')\V\xp(t')) 
-oo Jt J 



(21) 



x ( Xp (t) \Hi(t)\ * (<)> 



(22) 



It consists of three time-ordered steps by the electron: 
tunnel ionization, propagation in the laser field, and elas- 
tic scattering with the parent ion. Note that the SFA2 
used here is identical to the so-called improved strong 
field approximation 41, 42| . 

To evaluate the SFAl amplitude, we rewrite Eq. (QI 

as 



/i(k) = -i 



1 



dtE(t), 



iS(k,t) il p t 



(2tt) 3 / 2 J_ c 
J dre- i[k+A( * )l r rcos6»*o(r) 



(23) 



where I p is the ionization potential of the ground state 
^ / o( r ), and 6 is the polar angle. The ground state wave- 
function is calculated from the model potential V(r). To 
perform integration over space coordinates in (|23p . we 
use the identity 



(24) 



lm 



where ji (qr) is the spherical Bessel function and cos 8 = 
'Air I '3Yio(f) . Consequently, the integral over space co- 
ordinates can be expressed 

1<o (q) = J dre" iq r rcos6»* (r) 

= ^yyX^'yjmCS) J drr Z R nol() {r)]i{qr) 

lm 

x / dxY{ m {v)Y w {v)Y lomQ {v) (25) 



where the initial state wave function \&o( r ) = 
Rnoio( r )Yi m (r)- Due to reason mentioned before, for 
linearly polarized laser field, we consider mo = only, 
and 

J df^;„(f)y 10 (f)r ;omo (f) 

/ 3f9/„ 4- 1 > 

-C(lZ /;000)C , (lZ Z;000)5 m0 (26) 



3(2Z -I 


- 1 ) 


4tt(2Z - 


HI) 



where the C's are the Clebsch-Gordan coefficients. The 
remaining integration over r is done analytically if hy- 
drogenic wavefunction is used, otherwise it is evaluated 
numerically. The integration over time is carried out nu- 
merically. 

For the SFA2 amplitude, we used saddle point approx- 
imation for the integral over the momentum p of the 
intermediate states and Eq. (f22|) becomes 



/OO r-t 
dt / dt' 
-oo J — oo 



2tt 



-i 3/2 



e + i(t-t') 



E(t')e iIrt ' 



x e -i[S(P 3 ,t)-s(k,t)] e iS(p 3 ,t') 

X 



^-L^ J dre-^+^Kcos^oCr). (27) 

The saddle point is calculated with respect to quasiclas- 
sical actions only 

p s (t,t') = -j^ p ^dt"A(t"), (28) 

and the related actions are given by 

S(p s ,t) = l f Jt" [p^t') + A(t")f , (29) 
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and 



S(p s ,t') = l fjt" [p s (t,t')+A(t") 



(30) 



The arbitrary small parameter e in (|27p is introduced 
to remove possible singularity when t — > t' . Actually, the 
integral (|27|) converges for the case of initial state \&o( r ) 
having S symmetry (as for H) while it is divergent for 
the case of P symmetry (as for Ar) without e [43]. The 
Fourier transform of the potential V(r) in (9) is given by 



V{(\) = J dr exp(iq ■ r)V(r) 



-47T 



2a3ft4 



a 2 + q 2 



(31) 



It is obvious from (j3Tj> that the Fourier transform of V(r) 
diverges when q — > 0. Therefore, in actual numerical cal- 
culations, we multiply the potential by a damping factor 



V(r) = V(r)e~ ar 



(32) 



to avoid the singularity of the integral. We checked that 
the results mainly affect the magnitude but not the shape 
of the HATI spectra. 



D. Elastic differential elastic cross sections 

In this section, we briefly summarize the standard po- 
tential scattering theory which has been well documented 
in the textbooks, see Ref. [3|4^|, for example. Without 
loss of generality, here we address elastic scattering of 
an electron by a spherical potential V(r) by solving the 
time-independent Schrodinger equation 



[V 2 + k 2 - U(r)]^(v) = 



(33) 



where U(r) — 2V(r) is the reduced potential and k is the 
electron momentum, related to the incident electron en- 
ergy by k = \/2E. For short-range potential which tends 
to zero faster then r~ 2 as r — > oo, the scattering wave 
function satisfies the asymptotic outgoing wave boundary 
condition 



^+(r) r 



1 



(2tt) 3 / 2 



exp(ifcz) + f(0) 



exp(ifer) 



(34) 



where 9 is the polar angle measured from the incident 
direction. We choose the 2-axis along the direction of 
the incident wave vector k. 

We solve (13"3"|) by expanding the scattering wave func- 
tion into partial waves, 



ip+(r) 



2 1 

ir kr 



where Yi m is a spherical harmonic. The continuum waves 
are normalized to 8(k — k'). The radial function ui(kr) 
satisfies 



dr 2 



k 2 - 



1(1 + 1) 



-U(r) 



ui(kr) = 0. 



(36) 



For a plane wave, when U(r) — 0, the radial component 
ui(kr)/kr in (|35| is the standard spherical Bessel function 
ji(kr). 

The radial part of the scattering wave, ui(k, r), has the 
asymptotic form 



ui(kr) 



sin(A;r In 

v 2 



(37) 



where the phase shift Si reflects the influence of the in- 
teraction. 

The above equations are valid for short-range po- 
tentials only. For an electron in a Coulomb potential 
V c = —Z/r, its full wavefunction can be expanded as 



^c + (r) 



2 1 

7r kr 



Y,i l e la 'uUkr)Y lm (r)Y* n (k) (38) 



where 



ai = arg[r(; + 1 + it))] 



(39) 



is called the Coulomb phase shift with r] = —Z/k. The 
radial wavefunction u^(kr) is solved from 



dr 2 



1(1 + 1) 2rjk 



u^kr) = 0. 



(40) 



However, the expansion (|38|) does not converge well for a 
long-range Coulomb potential. For pure Coulomb scat- 
tering, the treatment in parabolic coordinates is simpler 
and the scattering amplitude is given by 



f c (0) = -rjexp(2ia 



cxp{-if7ln[sin 2 (#/2)]} 
1 2fcsin 2 (6>/2) ' 



(41) 



For electron-atomic target ion scattering within the 
single active electron model, the model potential is writ- 
ten as the sum of the Coulomb potential with Z = 1 and 
a short-range potential V s (r), see Eq. (jTUJ) ■ For such a 
modified Coulomb potential problem, the scattering am- 
plitude is given by 



f(o) = f c (e) + f(e) 



(42) 



where the first term is the scattering amplitude by the 
Coulomb potential alone [Eq. (|4"Tj) ] , and the second term 
is given by 



5^i' e «'t»,(*r)51m(*)^(k) (35) /» = E 



21 + 1 



exp(2ier ; ) exp(?<5;) sin<5 ; P/(cos#) (43) 



1=0 
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where the Pi(cos8) are Legendre polynomials, and Si is 
the phase shift from the short-range potential. Due to 
the short range nature, the summation in Eq. (|4"3"| can be 
truncated after some number of partial waves, depending 
the electron energy. The elastic scattering DCS is then 
given by 



a{k,6) 



dP 



|/c(0)+/(0)| 5 



(44) 



For high-energy collisions, one may calculate the dif- 
ferential cross sections using the first Born approxima- 
tion, or the plane wave Born approximation (P WB A) , in 
which, the DCS is given by 

ctpwba(M) = 4^(<l)| 2 ( 45 ) 

where q is the momentum transfer and its magnitude is 
q = 2k sin(0/2). In PWBA the continuum electron wave- 
functions are represented by plane waves. For electron- 
target ion collisions, PWBA is not valid even at large 
collision energies since it neglects the effect of long-range 
Coulomb interaction as well as the strong short-range po- 
tential due to the atomic ion or molecular ion core. In 
SFA2, scattering of the returning electron wave packet 
by the target ion is treated by the plane wave Born ap- 
proximation. This limits the accuracy of using SFA2 in 
describing the HATI spectra. 



III. QUANTITATIVE RESCATTERING THEORY 
AND ITS REGION OF VALIDITY 



A. Features of ATI electron energy and momentum 
spectra 

In the ATI photoelectron energy spectra, it has been 
well recognized that, after a sharp decrease to around 
3t/ p , a plateau exists from 4 to 10U p . This universal phe- 
nomenon has been observed in experiment, and in the 
TDSE and SFA calculations as well, see Fig. 2. These fig- 
ures also show that above about 4U P , SFA2 dominates the 
total electron spectra. Since SFA2 contains a first-order 
interaction between the electron and the target ion [see 
Eq. (I22p ], it is appropriate to attribute that ATI electrons 
above 4ZI p are produced by the rescattering processes. In 
Fig. 2 we note that the HATI spectra for Xe target ob- 
tained from TDSE are quite different from H when they 
are exposed to the same laser pulse. The HATI plateau 
in Xe is very flat. It remains almost constant between 
4.5-10C/p, while for H target, in the same energy region 
the yield drops by a large factor. 

Since ATI electrons are produced mostly along the di- 
rection of laser polarization (we consider linear polariza- 
tion only), at photoelectron energies where rescattering 
becomes dominant depend on the angle of the photoelcc- 
trons. In Fig. 3, we show the electron energy distribu- 
tions for electrons emitted at 0, 30, 60 and 90 degrees 
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FIG. 2: (Color online) Angle-integrated energy spectra (in 
units of the ponderomotive energy U p ) calculated from SFA 
(SFA1 and SFA2) compared with those by solving the TDSE 
for single ionization of (a) H and (b) Xe in a 5 fs laser pulse 
at the intensity of 1.0 x 10 14 W/cm 2 with the wavelength of 
800 nm. 
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FIG. 3: (Color online) SFA1 and SFA2 energy spectra for 
single ionization of Xe in a 10 fs laser pulse at the intensity 
of 1.0 x 10 14 W/cm 2 with the wavelength of 800 nm at angles 
of (a) 0°, (b) 30°, (c) 60°, and (d) 90°. The SFA2 dominant 
region moves to lower energies as the angle is increased. 



with respect to laser polarization, calculated using SFA1 
and SFA2, respectively. We note that the cutoff for SFA1 
where SFA2 becomes dominant shifts from about 4U P at 
zero degree to 3U P , 2U P , 1.5U P , respectively, at the angles 
of 30, 60 and 90 degrees. The laser parameters used in 
the calculations for Figs. 2 and 3 are given in the cap- 
tions. 

Another method of presenting the angular dependence 
of electron energy distribution is to display the 2D mo- 
mentum distributions. In Fig. 4, the 2D electron mo- 
mentum distributions obtained from TDSE and SFA2 are 
shown using the same laser parameters of Fig. 2. Only 
the large momentum portion is considered since the in- 
side is dominated by SFAl, i.e., the direct ionization. 
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photoelectron has momentum components 
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FIG. 4: (Color online) 2D electron momentum distributions 
for single ionization of H and Xe in a 5 fs laser pulse at the 
intensity of 1.0 x 10 14 W/cm 2 with the wavelength of 800 nm. 
(a) TDSE results for H; (b) TDSE results for Xe; (c) SFA2 
results for H; (d) SFA2 results for Xe. 



First we note that the appearance of circular "bands" 
at large momenta for the H target, in both calculations. 
The centers of these semi-circular rings are not at the 
origin, but are shifted along the polarization axis, one 
on each side. The rings are very similar for H and Xe 
in SFA2, but in TDSE, the intensity distributions show 
clear structure in Xe, in particular, clear minima at some 
angles. 



B. Extracting electron-target ion elastic differential 
cross sections from HATI spectra 



1993 



High-e nerg y ATI electrons have been observed since 



ii 



1 121 . 1131 ] . They were interpreted as 
due to the backscattering of the returning electrons by 
the target ion. Indeed classical simulation [l3j shows 
that electrons which return with maximum kinetic en- 
ergy of 3.17[/p, if backscattered by 180 degrees, would 
emerge with kinetic energy of about 10U p . Let Aq be 
the peak value of the vector potential of the laser pulse, 
Up = Aq/4. For an electron that returns at 3.17U P , it 
has momentum k r — \.26Aq. For a beam of electrons 
with momentum fc r , after elastically scattered, the mo- 
mentum space forms a circle in 2D (or a surface in 3D) of 
radius k r . Since scattering occurs when the laser field is 
nearly zero and the vector potential almost has the maxi- 
mum value Aq , each electron will gain an additional drift 
momentum Ao as it emerges from the laser field. These 
electrons were called back rescattered ridge (BRR) elec- 



trons in Morishita et al. [24(. The BRR electrons lie on 



a shifted circle in the photoelectron 2D momentum spec- 
trum. Let the direction of laser polarization be along 
the z-axis, and the y-axis perpendicular to it. After the 
backscattered electron emerges from the laser field, the 



k z = k cos 8 = ±Ao =F k r cos ( 
k y = k sin 8 = k r sin 9 r . 



(46) 
(47) 



The upper signs in Eq. (|46|) refer to the right-side 
(k z > 0) while the lower ones to the left-side (k z < 0). 
For backscattering, the angle 9 r is greater than 90°. 
These two equations can be expressed in vector form 
k = ±Aqz + k r . This vector relation (the outermost 
half-circle) is shown in Fig. 5 where the momentum is 
measured in units of Aq, and the angles are defined as 
shown. For k z > 0, the returning electron enters the tar- 
get from the right. After a large angle scattering, it is 
deflected by an angle 8 r . As the electron exits the laser 
field, it makes an angle 8 with respect to the polarization 
axis. 

In Morishita et al. [24] , it was argued that if the rescat- 
tering picture is correct, the intensity of electrons along 
the BRR should be proportional to the elastic DCS of 
the target ion by electrons with incident momentum k r . 
In 24|, this model was tested based on the HATI elec- 
tron spectra calculated from solving the TDSE, for H, 
Ne, Ar, and Xe targets. For the rare gas atoms, each 
target is represented by a model potential of the form, 
Eq. (9). The same model potential was used to calculate 
HATI electron momentum spectra D(k,8) and the elas- 
tic DCS, a(k r ,8 r ). By comparing the normalized yield 
of D(k, 8) along the ridge of fc r =1.26 A a with the DCS, 
a(k r ,8 r ), it was shown that the two indeed agree very 
well for the targets tested. The tests have been carried 
out for different laser intensities and mean wavelengths. 
Since the D(k,8) calculated from solving the TDSE are 
considered "exact", the test establishes the validity of 
attributing HATI electron momentum spectra to elastic 
backscattering of the returning electrons with momentum 
k r =1.26A . 



The theoretical result of [24{ has been limited to BRR 
electrons only where k r =1.26Ao, thus it leaves out a large 
portion of the HATI spectra where k r < 1.26Ao. For 
these lower energy electrons, as shown in Section II. A, 
the electrons may return to the ion core by following 
a long- or a short-trajectory. The ratio of the return- 
ing electron momentum k r (t r ) vs the vector potential 
A r =A(t r ) at the time of return, t — t r , are shown in 
Fig. 1(b). Note that the ratio does not change signifi- 
cantly in the time window for those returning electrons 
that can be backscattered to emerge with energies higher 
than 4Z7 p . Thus we set the relation k r = 1.26 A r for all 
HATI electrons, i.e., Eqs. ptj|) and (jl?)) are generalized 
to 

k z = kcosd = ±k r /1.26 =F K cos8 r , (48) 
k y = k sin 8 = k r sin 8 r . (49) 

Recall that if one neglects the effect of core potential, the 
returning electron momentum should be determined by 
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FIG. 5: (Color online) Typical 2D electron momentum distri- 
butions (in logarithmic scale): The TDSE calculation is for 
single ionization of Ar in a 5 fs laser pulse at the intensity of 
1.0 x 10 14 W/cm 2 with the wavelength of 800 nm. Photoelec- 
trons of a given energy are represented on a concentric circle 
centered at the origin. The elastic scattering of a returning 
electron with momentum k r in the laser field is represented 
by a partial circle with its center shifted from the origin by 
A r = fc r /1.26. High-energy plateau electrons are obtained via 
large-angle backscattering only. See text. 

the difference of the vector potentials at the return time 
and the born time. From Fig. 1(a), the range of born time 
is very narrow, thus the same relation between k r and A r 
for all HATI electrons is a good approximation. With this 
model, the center of the circle for each momentum k r is 
shifted by /c r /1.26. This has important implications since 
the laser parameters such as peak intensity or wavelength 
do not enter explicitly in Eq. (|48|) or (|49|) any more. How 
good is this model? We test its validity using accurate 
numerical results from TDSE calculations. 

Recall that the validity of Eq. flU, Eqs. gSJ) and (|39). 
has been fully tested for the case of short pulses and for 
k r = 1.26^0 m 01 ■ I n Fig- 6( a )7 we show the theoret- 
ically calculated DCS for Ar at k r — 1.22 and compare 
with the DCS extracted from HATI momentum spectra 
obtained from solving the TDSE for Ar target. In one 
case we use a five-cycle pulse with 800 nm mean wave- 
length, and peak intensity of 1.0 x 10 14 W/cm 2 . This is 
the same as the BRR discussed in In another case, 
the peak intensity used is 1.4 x 10 14 W/cm 2 . The BRR 
electron momentum for the latter is 1.32. If we extract 
the DCS from the HATI spectra of the latter at fc r =1.22, 
as seen from Fig. 6(a), the results are still quite good. 

Up to now, we have focused on short laser pulses. For 
longer pulses, the electrons generated from different op- 
tical cycles of the laser can interfere. Such interference 
would result in the well-known ATI peaks which are sep- 
arated by the photon energy of the laser. From Fig. 5, 
along £; r =constant, the photoelectron energy changes as 
9 r is varied. Thus the assumption that W(k r ) is con- 
stant along a fixed k r is no longer correct because of the 
interference in the wave packet. However, the interfer- 
ence effect is well behaved. In Fig. 6(b), we show the 
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FIG. 6: (Color online) Electron yield extracted for fixed k r 
from the 2D electron momentum distributions for Ar and Xe 
by solving the TDSE compared with the corresponding elastic 
scattering DCS's (thick black solid lines), (a) For the target 
of Ar and k r = 1.22 by a 5-cycle laser pulse with the wave- 
length of 800 nm at intensities of 1.0 and 1.4 x 10 14 W/cm 2 , 
respectively; (b) For the target of Xe and k r — 0.92 by a 8- 
cycle laser pulse with the wavelength of 800 nm at intensities 
of 5.8 , 6.2 and 7.0 x 10 13 W/cm 2 , respectively. 

DCS for e-Xe+ at fc r =0.92. We also extracted the DCS 
at the same A; r =0.92 using the HATI spectra generated 
by three lasers, with peak intensity of 7.0, 6.2, and 5.8 in 
units of 1.0 x 10 13 W/cm 2 , for 800 nm lasers of durations 
of 8 cycles. By assuming a constant wave packet W(k r ), 
the "extracted" DCS oscillates but the peak positions of 
the oscillation still follow the DCS quite accurately (after 
normalized) . Such oscillations appear to be worrisome in 
attempts to extract DCS from HATI spectra using longer 
laser pulses. However, as will be shown in Section IV. D, 
this is not a problem for extracting DCS from exper- 
imental HATI spectra since experiments "intrinsically" 
integrate electron spectra generated over a distribution 
of laser intensities. More discussions on the oscillations 
of the wave packet are given in the next subsection. 

C. Extracting rescattering wave packet from the 
HATI spectra 

Given the relation between (fc, 9) and (k ri 9 r ) in 
Eqs. (|48|) and (|49|. in general one can write D(k,9) = 
W(k r ,9 r )u(k r ,9 r ). If the rescattering concept is mean- 
ingful for a given k r , we expect W(k r ,0 r ) = W(k r ), i.e., 
the rescattering wave packet distribution is independent 
of scattering angles. To illustrate this point, we calculate 

W(k r ,9 r ) = D(k,9)/a(k r ,6 r ), (50) 

where D(k, 9) is obtained from TDSE and <r(k r , 9 r ) from 
Eq. (|44p. The results for W(k r ,9 r ) are shown for the 
angular range of 9 = 155° to 180° in Figs. 7(a) and 
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FIG. 7: (Color online) Right-side wave packets(fc z > 0) ex- 
tracted from the TDSE and SFA2 electron momentum distri- 
butions for single ionization of Ar in a 5 fs laser pulse with 
the wavelength of 800 nm. (a,c) TDSE and SFA2 results at 
the intensity of 1.0 x 10 14 W/cm 2 . (b,d) TDSE and SFA2 
results at the intensity of 2.0 x 10 14 W/cm 2 . 



7(b), for two intensities (with other laser parameters 
given in the figure), respectively. It is quite clear that 
there is little angular dependence of 9 r such that we can 
write W(k r ,6 r ) = W(k r ). We emphasize that this re- 
lation is based on computational results where D(k,0) 
and a(k r ,9 r ) are calculated "exactly". The correctness 
of W(k r , 9 r ) = W{k r ) justifies the relation in Eqs. 
and (|49]1 and it provides a strong statement of the QRS 
model for HATI electron momentum spectra, as stated 
in Eq. (Q]) in the Introduction. Note that W(k r ) is ex- 
tracted from D(k, 0) at the end of the laser pulse, thus 
W{k r ) includes all the quantum interference due to the 
long- and short-trajectory electrons, and interference due 
to wave packets generated from different optical cycles. 
In Figs. 7(a,b), the pulse duration is 5 fs. The oscillation 
in the wave packet is due to interference of long- and 
short-trajectory electrons that return with the same k r . 
As the laser intensity increases, the oscillations become 
faster. This increase of oscillations can be easily under- 
stood based on the SFA2. Note that there are two wave 
packets, W(k r ), one from the left and the other from the 
right toward the target. 

As discussed earlier, in strong field approximation, 
SFA2 dominates over SFA1 for HATI electrons with en- 
ergies above 4£/ p , see Fig. 2. Since rescattering is in- 
cluded in SFA2, we also check whether separation simi- 
lar to Eq. |T]) also applicable to D(k, 8) calculated from 
SFA2. In SFA2 the elastic scattering of RWP with the 
target ion is treated to first order only, thus the cor- 
responding a(k r ,8 r ) is calculated using the plane-wave 
Born approximation, Eq. (|45[) . Following the same pro- 
cedure as for the TDSE results, we extract the RWP from 
SFA2. The results are shown in Figs. 7(c) and 7(d) for 
the same two sets of lasers used in the TDSE calcula- 
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FIG. 8: (Color online) (a) Right-side wave packets (k z > 
0) extracted from TDSE for single ionization of Ar in a 5- 
cycle pulse with the wavelength of 800 nm at the intensity of 
l.OxlO 14 W/cm 2 at r = 16O°, 170° and 180°, respectively; (b) 
Same as (a) but for an 8-cycle pulse; (c,d) Electric field and 
vector potential used in (a,b), respectively. The peak values 
of both E(t) and A(t) are normalized to 1. The horizontal 
lines in (c,d) are for A r =0. 625^4.0. 



tions. Comparing Figs. 7(a,c) and Figs. 7(b,d), respec- 
tively, we note that the RWP's extracted from SFA2 and 
from TDSE have very similar shape. After normalization, 
their dependence on k r is nearly identical for the same 
laser. The absolute value of W(k r ) from SFA2 is smaller 
in general since ionization yield calculated using strong 
field approximation is smaller in general. The similarity 
of W(k r ) from SFA2 and from TDSE also allows us to 
interpret the increase of oscillations of RWP as the laser 
intensity is increased, as seen in Fig. 7. This increase 
can be traced to the actions S, Eqs. ([29]) and (f30|) . whose 
values increase quadratically with the vector potential, 
or linearly with the ponderomotive energy U p and time 
duration of the pulse. 

For longer pulses, the RWP should reflect the inter- 
ference from different optical cycles. Consider the wave 
packets shown in Figs. 8(a,b), generated by an 800 nm, 
peak intensity of 1.0 x 10 14 W/cm 2 , but one with 5, and 
the other 8 optical cycles, respectively. The correspond- 
ing .E-field and A- vector are shown to the right. The 
horizontal lines in Figs. 8(c,d) are for A r = 0.625^4o- 
Electrons return with this A r have k r = 1.26 A r that 
can be back rescattered to reach HATI energy of 4?7 p . 
Note that A r determines the photoelectron energy, but 
the yield is determined by the electric field at about 3/4 
cycles earlier. Thus for the 5-cycle pulse, HATI electrons 
are generated from one optical cycle only, and the oscil- 
lation seen in Fig. 8(a) is due to interference from long- 
and short-trajectory electrons. For the 8-cycle pulse, at 
least two optical cycles make contributions to the HATI 
spectra, thus interference seen in Fig. 8(b) becomes much 
more numerous. The fast oscillations from electrons gen- 
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erated at different optical cycles overwhelm the slower 
oscillations from electrons generated within the same cy- 
cle, such that the long- and short-trajectory interference 
is seen as the oscillation of the envelope in W(k r ). As 
the pulse duration increases, the momentum distribution 
of the RWP will become flatter except for the cutoff re- 
gion. As demonstrated in Fig. 8(b), the finer oscillations 
in the RWP extracted from different angles 9 r will not be 
the same. However, the envelope of the wave packet is 
independent of 6 r , such that a single wave pac ket is still 
meaningful. More examples can be seen in 



D. Target independence of the shape of the 
returning electron wave packet 



The fact that the W(k r ) extracted from SFA2 is similar 
in shape to that extracted from TDSE is again a conse- 
quence of the validity of the rescattering model. While 
the absolute returning electron yield is determined by the 
initial tunnel ionization rate, its momentum distribution 
W(k r ) is determined almost entirely by the laser field. 
In SFA2, this interaction is fully included. The electron- 
target ion interaction, which is included in TDSE, affects 
W{k r ) weakly only since the returning electron spends 
most of the time away from the target ion where the 
field is dominated by the laser's electric field. 

The fact that W(k r ) can be obtained from SFA2 rea- 
sonably accurately has a far-reaching implication. Before 
discussing such implications, first we look more carefully 
at how well the W{k r ) obtained from TDSE and SFA2 
agree, for different target atoms. In Figs. 9(a-d) we show 
these RWP's generated by an 800nm, 8-cycle, peak in- 
tensity 1.0 x 10 14 W/cm 2 laser pulse for H, Ar and Xe 
targets. We show the "left" and "right" RWP's which 
are different for the short pulse used. To first order, all 
the wave packets on the left are similar, and all the wave 
packets on the right are similar. A more careful exam- 
ination reveals that there are differences. Among the 
different targets calculated using TDSE, we note that 
the distribution of W(k r ) tends to shift to higher k r as 
the ionization energy increases. This is expected since 
the returning electron is seeing a more attractive poten- 
tial from the ion core, and the effect is bigger for lower 
energy electrons. Still the effect is only a few percents. 
For longer pulses or higher intensities, see Figs. 9(e,f), 
the wave packets at lower momenta are mostly very flat 
where small shift of individual peaks is not very impor- 
tant. Due to the small difference of the RWP on the 
target for a given laser pulse, the wave packet can be 
generated using SFA2 from a hydrogenlike target, with 
the effective charge chosen so that it gives identical bind- 
ing energy of the target. 
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FIG. 9: (Color online) (a) Left-side wave packets ( k z < 0) 
extracted from TDSE for single ionization of H, Ar and Xe 
in a 8-cycle pulse at the peak intensity of 1.0 x 10 14 W/cm 2 
with wavelength of 800 nm, respectively; (b) Same as (a) but 
for the right-side (k z > 0); (c) Same as (a) but from SFA2; 
(d) Same as (b) but from SFA2; (e) Right-side wave packets 
extracted from SFA2 for single ionization of H and Ar in a 
5-cycle laser pulse at the intensity of 1.0 x 10 14 W/cm 2 with 
the wavelength of 800 nm; (f) same as (e) but for intensity of 
2.1 x 10 14 W/cm 2 and wavelength of 2000 nm. 



E. Elastic electron-ion differential cross sections at 
large angles 

Elastic and inelastic scattering cross sections between 
free electrons and atomic ions have been studied in cross- 
beam or merged-beam experiments [46i j . In such exper- 
iments, a well-collimated electron beam with precisely 
defined energy is prepared. For neutral atomic or molec- 
ular targets which can be placed in a gas cell, there have 
been lots of experimental and theoretical investigations 
in the past half a century. These studies tend to focus 
on sharp features like Feshbach resonances which often 
require a careful treatment of electron correlation effects. 
For the HATI spectra, the returning electron is described 
by a wave packet which has a broad momentum distri- 
bution. Here we tend to focus on the broader energy and 
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FIG. 10: (Color online) (a) DCS's for Xe at k r = 0.8 in the 
angular range of 110°-180°. Solid line: "exact" scattering 
wave result [Eq. (42)]; broken line: PWBA [Eq. (45)]; Chain 
line: Coulomb wave result [Eq. (41)]. (b,c,d) 2-dimensional 
"exact" scattering wave results of the DCS's for Ar, Kr and 
Xe, respectively. 

the angular dependence, and neglect the many-electron 
effects to first order. Within this model, the calculation 
of elastic scattering cross sections is quite simple, as de- 
scribed in Section II. D. 

In Fig. 10(a) we compare the DCS's calculated for Xe 
at the incident momentum of k r = 0.8 for scattering an- 
gles from 110°-180°, where the continuum wave functions 
are represented by "exact" scattering waves, by plane 
waves and by Coulomb waves, respectively. The DCS's 
are normalized near about 130°. We note that the DCS 
calculated from the PWBA is rather featureless, so is 
the DCS calculated using Coulomb wavefunctions, each 
drops monotonically with increasing scattering angle. On 
the other hand, the DCS calculated from the scattering 
wave shows complicated pattern which is the well-known 
Ramsauer-Townsend electron diffraction. This illustrates 
that both plane wave and Coulomb wave are very poor 
approximation for describing electron-ion collisions, espe- 
cially for electrons which undergo large-angle scattering. 
For such large deflection angles, the electron has to pen- 
etrate the ion core, thus seeing the short-range part of 
the potential. Without scattering waves, the diffraction 
by the strong potential V s (r) is neglected and the HATI 
spectra cannot be correctly reproduced. 

In Figs. 10(b-d) we show a(k r ,9 r ) for Ar, Kr and Xe 
for k r from 0.6 to 1.1. The complicated structure, and 
the increase of cross sections close to 180°, are quite ev- 
ident for all three targets. Such structure would not ap- 
pear if H is used as the target. For the large angles 
discussed here, the minima in the DCS come from the 
interference of contributions from several partial waves 
in the scattering by the short range potential. At large 
angles, the Coulomb scattering amplitude [See Eq. (|4ip] 



is relatively small. The interference of the two ampli- 
tudes in Eq. (I42|) can also produce interference minimum 
shown in typical textbook examples [3], but such min- 
imum occurs at smaller angles. (See the two minima in 
the DCS for e~+Ar+ collisions in Fig. 4(c) of [28| where 
the large-angle minimum at 140° was derived from the 
HATI spectra and the small-angle minimum at 89° was 
observed from electron-Ar + colliding beam experiment 
0.) 

The strong dependence of a(k r ,6 r ) on the target po- 
tential shows that HATI spectra can be obtained accu- 
rately only if electron-ion scattering is treated accurately 
in the nonlinear interactions of lasers with atoms and 
molecules. Thus TDSE calculations using approxima- 
tions where the singularity of the Coulomb core potential 
is regularized should be handled with care. Such regular- 
ized models have less effect on the total ionization yield 
or the electron spectra at low energies, but it will affect 
the HATI spectra since these electrons undergo close col- 
lisions with the target core. This also implies that strong 
field calculations where continuum electrons are treated 
at the level of plane waves (as in SFA2), Coulomb waves 
(such as using Coulomb Volkov states), or in the eikonal 
approximation, will not be adequate for the description 
of the HATI spectra. 

F. The practical quantitative rescattering model and 
its validity 

Based on the established validity of the QRS model 
using TDSE and SFA2, we now propose a practical QRS 
model for obtaining HATI electron momentum spectra. 
We would apply this model to electron energies above 
about 4Up for the total electron energy spectra. For the 
energy dependence at each fixed angle, the lower energy 
limit where this theory applies can be relaxed, as illus- 
trated in Fig. 3. 

In the practical QRS model, the HATI momentum dis- 
tribution D(k,8) is calculated using Eq. ([TJ) . We obtain 
the DCS using Eq. (|44]) and the returning electron wave 
packet from 

W(k r ) = D SFA2 (k, 9)/a PWBA (k r ,e r ), (51) 

where Z)sfA2(^,#) = I/2I 2 and f 2 is calculated from 
Eq. (f27|) . and <JpwBA(k r , r ) is calculated from Eq. (|45j) . 
The relation between (k,9) and (k r ,8 r ) are given by 
Eqs. (gHl) and (|49]). To obtain W(k r ) from Eq. (HO}, we 
need to do SFA2 calculation at one angle 9 r only. We 
typically use 6> r =170°. 

Since W(k r ) is nearly independent of the target for a 
given laser pulse, we can also perform SFA2 calculation 
using hydrogenic potential with an effective charge that 
reproduces the binding energy of the target atom. We 
emphasize that using the QRS, the absolute yield is not 
obtained. This is also true for most of the energy and 
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FIG. 11: (Color online) Logarithmic photoelectron 2D mo- 
mentum distributions by a 5-cycle pulse at the peak intensity 
of 1.0 x 10 14 W/cm 2 with wavelength of 800 nm. (a) TDSE 
results for Ar; (b) QRS results for Ar; (c) TDSE results for 
Xe; (d) QRS results for Xe. 



momentum spectra reported in experiments. Using the 
practical QRS model, the calculation of HTAI spectra 
can be a factor of thousands faster compared to TDSE 
calculations. However, this is useful only if we can show 
that the QRS reproduces the HATI momentum and en- 
ergy spectra at the level comparable to TDSE results. 

In Fig. 11 we show the 2D electron momentum spectra 
calculated using QRS and TDSE for Ar and Xe using the 
laser parameters indicated in the figure. Recall that for 
short pulses the left and right wave packets have to be 
calculated separately. One can see the overall agreement 
between the QRS and TDSE calculations. In making 
color plots we renormalize the spectra and the same color 
schemes are used in the figures. Since the same laser is 
used in the calculations for both Ar and Xe, according 
to the QRS, the difference in Figs. ll(c,d) are almost en- 
tirely due to the difference in the DCS [see Figs. 10(b,d)] 
by the returning electrons colliding with the target ions 
of Ar and Xe, respectively. 

To display the comparison quantitatively, we show the 
total electron energy spectra obtained from TDSE and 
those from QRS for energies above AU p . In Fig. 12(a), 
the electron energy spectra are obtained by integrating 
the momentum spectra in Fig. 11. The QRS and TDSE 
results agree quite well above AU p . For Ar, we see some 
discrepancy close to AU P . In Fig. 12(b), we compare the 
total electron spectra using lasers of three different wave- 
lengths, with the same number of cycles, but with the 
intensity adjusted such that U p — 6 eV for each pulse, 
and the Keldysh parameter is 1.14. The electron spec- 
tra calculated from TDSE arc placed on the absolute 
scale, while the QRS results arc normalized individually 
to achieve best agreement. Among the three cases, the 
DCS's are the same. Even though the returning electron 
wave packets cover the same range of momentum, the 



FIG. 12: (Color online) Comparison of angle-integrated pho- 
toelectron energy spectra from TDSE and QRS. The TDSE 
spectra are plotted from to 12U P and the QRS results are 
shown from 4 to 12U P . The QRS spectra are relative and 
normalized to TDSE results at high energies, (a) Single ion- 
ization of Ar and Xe in a 5-cycle pulse at the peak intensity 
of 1.0 x 10 14 W/cm 2 with wavelength of 800 nm; (b) Single 
ionization of Ar in a 5-cycle pulse at peak intensities of 4.0, 
1.0 and 0.64 x 10 14 W/cm 2 with wavelengths of 400, 800 and 
1000 nm, respectively. 



interference features in W(k r ) are different. Such inter- 
ference features are easily reproduced in the wave packet 
calculated from SFA2. For 400 nm, the TDSE results 
show clear ATI peaks which are not reproduced in the 
QRS model, as explained in Section III.D. For additional 
examples of the comparison between the QRS model and 
TDSE, see Chen et al. (25 1. 



IV. COMPARING QRS MODEL WITH 
EXPERIMENTAL ATI ELECTRON SPECTRA 



A. Volume effect in experimental ATI spectra 

In the previous Section we compared the HATI spectra 
calculated using TDSE and QRS for a given laser pulse 
with single intensity. Experimentally, the intensity dis- 
tribution of a focused laser beam is not uniform in space. 
The HATI electrons are collected from the whole focused 
volume. Thus to compare with experimental HATI spec- 
tra, theoretical calculations must include the volume ef- 
fects For a peak intensity Iq at the focal point, the 
yield of the photoelectrons with momentum k should be 



S(k,I )=P 



Ia fdV 



dl 



(52) 



where p is the density of atoms in the chamber, Di(k, 9) 
denotes the momentum distribution for a single inten- 
sity / and (dV/dI)dI represents the volume of an isoin- 
tensity shell between / and I = I + dl defined in 47 1 
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for a Lorentzian (propagation direction) and a Gaussian 
(transverse direction) beam profile. We use the trape- 
zoidal rule for the integration over intensity with suffi- 
ciently small step size of 0.01 x 10 14 W/cm 2 . 

In the QRS calculations, we obtain the volume- 
integrated returning electron wave packet using Eq. (Q]) 
since the DCS does not depend on the laser intensity. 
Consequently Eq. (|52|) becomes 

S(k,I ) = W Io (k r )a(k r ,6 r ) (53) 

where Wi (k r ) is the volume-integrated wave packet at 
the peak intensity Iq 

Wi (k r ) = pj*° (jjp) dl (54) 

with Wi(k r ) being the wave packet for the laser pulse at 
a single intensity I. 

B. Wavelength dependence 

In a recent experiment, Colosimo et al. [48| reported 
electron spectra generated from Ar by infrared (IR) to 
mid-infrared (MIR) lasers with the same peak intensity, 
for wavelength of 800, 1300, 2000 and 3600 nm, respec- 
tively. The measured electron spectra with electron en- 
ergies in units of U p are shown in Fig. 14, as well as the 
results from the QRS model above AU P where volume in- 
tegration effect has been included. First we note that 
there is a general agreement between the HATI spec- 
tra from the measurement and from the QRS model. 
For 800 nm, the QRS underestimates the 4-5.5C/ p region 
which could be due to the resonantlike enhancement ef- 
fect [49|]. In the figure, the electron spectra from the 
different wavelengths are normalized near at threshold. 
One notes that the HATI yields decrease rapidly with 
increasing wavelength. Such decrease is also familiar in 
the HHG spectra [50j, empirically estimated to decrease 
like A -5 5 . For longer wavelength lasers, the excursion 
distance and the return time of the electron after tunnel 
ionization both scale with the wavelength, thus increase 
the effect of broadening in the returning wave packet. 

From Fig. 13, we note that the slope of the electron 
spectra appears to flatten out considerably in going from 
1300 nm to 2000 nm. To understand the origin of the 
slope change, we show in Fig. 14(a) the returning elec- 
tron wave packet obtained from the QRS model. The 
momentum of the wave packet is expressed in units of Aq 
of the laser pulse at the focus center. We note the rapid 
oscillations are due to the ATI peaks. In Fig. 14(b) the 
volume integrated wave packets are shown. We note that 
they are mostly flat above 4[/ p . According to the QRS 
model, Eq. (53), we thus expect that the slope change 
in the experimental data is due to the elastic scattering 
cross sections. In Fig. 15, we show the differential cross 
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FIG. 13: (Color online) Comparison of experimental and QRS 
angle-integrated photoelectron energy spectra for single ion- 
ization of Ar by few-cycle laser pulses at the peak intensity 
of 0.8 x 10 14 W/cm 2 with wavelengths of 800 nm, 1300 nm, 
2000 nm and 3600 nm, respectively. The experimental mea- 
surements are taken from [43] • The QRS results starting from 
4U P are normalized to the experimental data individually at 
high energies to get best fit. Electron energies are expressed 
in units of U p . 
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FIG. 14: (Color online) Wave packets at 800 nm, 1300 nm, 
2000 nm and 3600 nm. (a) For single intensity of 0.8 x 10 14 
W/cm 2 ; (b) Volume-integrated for the peak intensity of 0.8 x 
10 14 W/cm 2 at the focal point of a Gaussian beam. See text. 



sections of electron-Ar + collisions in the four momen- 
tum ranges for the returning electrons that contribute to 
the HATI spectra for the four wavelengths used. Note 
that the DCS peaks sharply at large angles near 180° for 
the momentum range of 1.60-2.65 for the 2000 nm laser 
pulse. Such behavior of the DCS is responsible for the 
much flatter HATI spectra seen in Fig. 15 for the 2000 
nm pulse. Note the similarity of the DCS in Fig. 15(c) 
in Ar, and the DCS in Xe in Fig. 12(a). The large DCS 
at angles close to 180° is responsible for the much flatter 
energy dependence observed in the HATI spectra. 

The above analysis shows that the QRS model not only 
can explain experimental results, but also offers a clear 
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FIG. 15: (Color online) DCS's for elastic scattering of elec- 
tron with Ar + at large scattering angles of 110°-180°. The 
incident electron momenta are: (a) 0.65-1.10, (b) 1.00-1.75, 
(c) 1.60-2.65, and (d) 2.95-4.80, corresponding to the high 
energy plateau regions for 800 nm, 1300 nm, 2000 nm and 
3600 nm, respectively. 
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FIG. 16: (Color online) (a) Comparison of the total energy 
spectra of QRS with experimental measurements for Ar of flOl ] 
at different intensities; (b) Illustration of how to adjust inten- 
sity in the QRS simulation to fit experimental measurement. 
In the QRS calculations, volume effect has been included. 



physical interpretation of the origin in the difference of 
the observed HATI spectra. We comment that in princi- 
ple these HATI spectra can also be calculated by solving 
the TDSE. However, as shown in Colosimo et al. [Zjj ]. 
such TDSE calculations for MIR lasers are very difficult 
due to their large long excursion distance and the need 
of a large box in the calculation for MIR laser pulses. 
Using the QRS, the calculation is much easier. In fact, 
the volume integrated returning wave packet is quite flat 
that one may even just approximate it by a constant. 
The QRS model then would predict that the slope seen 
in Fig. 13 is due to the integration of the DCS over the 
scattering angles. Note that by changing the peak in- 
tensity at the laser focus, the same flat plateau seen in 
Fig. 13 is expected to shift to different wavelength. 
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FIG. 17: (Color online) Comparison of the total electron en- 
er gy s pectra from QRS with the experimental measurements 
of ( a ) Xe; (b) Kr. In the QRS calculations, volume effect 
has been included. 



C. Target and intensity dependence 

We also used the QRS model to simulate HATI spec- 
tra from some earlier experiments of Grasbon et al. 13] • 
Since in experiments the peak laser intensity is often de- 
termined approximately only, we perform the simulation 
by treating the peak laser intensity as a free parameter. 
We assume that all the electrons from the focal volume 
are collected in the experiment. We use the pulse du- 
ration and wavelength reported in the experiment. In 
Fig. 16(a) we compare the experimental HATI spectra 
measured with the ones obtained from the QRS model. 
Note that the peak intensity reported in the experiment 
is often different from the one that gives the best fit ob- 
tained from the QRS. In each spectra, we normalized the 
QRS result such that it gives best overall fit to the exper- 



imental data. Fig. 16(b) illustrates the example how the 
slope of the HATI spectra changes as the peak intensity 
at the focus is varied. 

Similar analysis has been carried out for Xe and Kr 
targets, with the results shown in Fig. 17. In Kr, the 
agreement is very good for all three spectra. The small 
difference at the low energy end could be due to contribu- 
tions from direct tunneling ionization. For Xe, we have 
been able to reproduce the outer part of the plateau well. 
At lower electron energies the experimental data show a 
deeper minimum for the two upper intensities which are 
not reproduced by the QRS simulation. The reason for 
this discrepancy is not clear at this time. In Fig. 18(a) 
we show the relative weights of the different peak inten- 
sities that contribute to the total ionization yields for the 
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FIG. 18: (Color online) Single ionization of Xe by a 8 fs pulse 
at the peak intensity of 0.65 x 10 14 W/cm 2 with wavelength of 
800 nm. (a) Laser volume effect analysis showing the relative 
contributions of laser peak intensity to the generated electron 
spectra, (b) Comparison of QRS with TDSE at the single 



Comparison of QRS with 



intensity of 0.59 x 10 i4 W/cm' 
TDSE with volume effect included. 



focal peak intensity of 0.65 x 10 14 W/cm 2 . In Fig. 18(b), 
we show that the HATI spectra from the QRS agree well 
with the result from solving the TDSE at peak intensity 
of 0.59 x 10 14 W/cm 2 , the intensity that contributes most 
to the electron yields, see Fig. 18(a). 

Fig. 18(c) shows that volume integrated electron spec- 
tra calculated from TDSE and from the QRS agree well, 
but both cannot reproduce the deeper experimental min- 
imum. This minimum occurs at about 20 eV, which is 
close to returning electron momentum of k r = 0.67, or 
electron energy of about 6 eV. At such a low energy, 
electron- Xe + elastic scattering cross sections may not 
be well described by the model potential approach. Re- 
call that in QRS and TDSE, we treated Xe atom using 
a single active electron model. For electron-atom and 
electron-ion collisions, it is generally known that many- 
electron correlation effect becomes more important as the 
collision energy decreases. (An example where the sin- 
gle electron approximation fails has been noted in the 
photodetachment of H~ by MIR lasers, see Zhou et al 
|43(.) Thus one possible explanation for the failure of 
the QRS theory to reproduce the experimental observa- 
tion in Fig. 17(a) is the need to include electron correla- 
tion effect in calculating e~-Xe + elastic scattering cross 
sections. While such calculations have been carried out 
using many-body perturbation theory 5^, 53[ in a num- 
ber of cases, no such results have been reported for the 
present system. The fact that the minimum occurs at 
the same photoelectron energy in the experimental data 
at the two upper intensities, see Fig. 17(a), also offer a 
hint that the discrepancy is due to error in the DCS used. 
The minimum does not appear at the lowest intensity in 
Fig. 17(a) since the HATI yield drops rapidly in the same 
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FIG. 19: (Color online) Comparison of the DCS's with elec- 
tron yield extracted from the experimental measurements for 
Xe by laser pulses with different CEP's and unlocked CEP 
for (a) fc r =0.9, (b) *v=1.0, (c) fc r =l.l, and (d) k T =1.2. The 
CEP's are given in the figure for phase- locked pulses. 



energy region. 



D. Extracting DCS from experimental electron 
spectra and Other Applications 

We have applied the QRS model to a number of other 
topics involving HATI electrons so far. In Chen et al. 25 1 
the flatness of the HATI spectra vs electron energies ob- 
served for rare gas atoms and alkali atoms were inves- 
tigated and interpreted in terms of the energy and an- 
gular dependence of the DCS, similar to the examples 
presented here. The QRS model has also been applied to 
retrieve the absolute value of the carrier-envelope-phase 
of few-cycle pulses, as well as the pulse duration and the 
peak laser intensity in Micheau et al. 0, [3(3]. Using 
the wave packet extracted from the HATI spectra and 
electron-impact ionization cross sections, the nonsequen- 
tial double ionization of Ar has been obtained in Micheau 
et al. 



31] 



Following the earlier theoretical paper of Morishita et 
al. [23] in which DCS was extracted from the HATI spec- 
tra along the BRR, i.e., along k r = 1.26 An, the prediction 
was confirmed in two experiments 2^, 27 1. In these ex- 
periments, pulse durations of about 100 fs and 8 fs were 
used, respectively. As shown in Fig. 6(b), for such long 
pulses the extraction of DCS is difficult because of the in- 
terference in the wave packet. Interestingly, since the ex- 
perimental electron spectra are collected from the whole 
focal volume, the oscillation in the volume-integrated 
wave packet at a given k r is smoothed out. Thus in ef- 
fect, the success of retrieving the DCS in these two earlier 
experiments is based on the present QRS model, 

rather than the early theory of [2^]. In other words, the 
effect of volume integration in the electron spectra ac- 
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tually simplifies the retrieval of the DCS. Similarly, the 
DCS can also be extracted from ATI spectra generated by 
few cycle pulses. For these pulses, the returning electron 
wave packet depends on the carrier-envelope phase. One 
can obtain the DCS from experiments where the CEP 
is locked, or from measurements where the CEP is not 
locked. The extracted DCS should be the same based 
on the QRS model. In Fig. 19, we show the extracted 
DCS for Xe at four different electron momentum values, 
using the experimental data of Kling et al. [51]. The 
data extracted from these different data sets agree quite 
well, and they agree reasonably well with the theoretical 
DCS calculated from single active electron approxima- 
tion. The theory shows slightly deeper minimum, but 
the angular resolution was not considered in the theoret- 
ical calculations. We emphasize that the extracted DCS 
should be independent of the lasers used. 



SUMMARY AND OUTLOOK 

In this paper we present a comprehensive QRS theory 
for describing the energy and momentum distributions of 
HATI electron spectra generated by intense laser pulses. 
Although HATI spectra have been interpreted in terms of 
rescattering concept since the 1990's, the QRS model is 
the first quantitative rescattering theory that can achieve 
accuracy comparable to those obtained from solving the 
TDSE. The essential ingredient of the QRS is governed by 
Eq. (1) which states that HATI electron momentum dis- 
tributions can be expressed as the product of a returning 
electron wave packet W(k r ) and the elastic differential 
cross sections a(k r ,8 r ) between free electrons and the 
target ion. The validity of the QRS model is carefully 
tested against results obtained from solving TDSE using 
atomic targets in the single electron approximation. Here 
are a number of the most notable results: 

(i) The wave packet W(k r ) is mostly independent of 
the target, and can be calculated from the 2nd-order 
strong field approximation theory. All the laser depen- 
dence in the HATI spectra is governed by W(k r ). The 
target structure enters through the DCS, <r(k r , r ). 

(ii) Using a practical QRS model, HATI spectra can be 
calculated with accuracy comparable to those obtained 
from solving TDSE, but with saving of computing time 
by a factor of thousands. This is particularly important 
when one wants to compare calculations with experimen- 
tal measurements where repeated calculations have to be 
done for many intensities. With the QRS, as we demon- 
strated in Section IV, quantitative comparison of HATI 
spectra between theory and experiment is now possible. 
The QRS theory also is capable of checking consistency 
of different experiments. 

(iii) The QRS model allows the extraction of the re- 
turning electron wave packet (in the laser field) based on 
HATI spectra at the end of the laser pulse, instead of us- 



ing classical simulations or theoretical models when the 
electron is still in the laser field. Our extracted W(k r ) 
reflects the interference due to returning electrons with 
long- vs short-trajectories, as well as interference from 
electrons released from different optical cycles. 

(iv) The QRS model allows one to separate the role of 
laser pulses which is contained in W(k r ), and the role of 
the target which is contained in a(k r ,6 r ). By calculat- 
ing W(k r ) using SFA2, the difficult part of the nonlinear 
laser-atom interaction is avoided. This accounts for the 
major saving in computer time in the QRS as compared 
to TDSE. The conceptual separation of the nonlinear ef- 
fect of lasers on the wave packet and the structural in- 
formation contained in the DCS has many ramifications. 
In particular, it should be possible to generalize the QRS 
model to include many-electron effects, just by replacing 
a(k r , 6 r ) with the DCS calculated including electron cor- 
relation effect. But most importantly, the QRS model 
offers an opportunity to calculate accurate HATI spectra 
for molecules. 

For molecular targets, there are many possible impor- 
tant applications. One can use a pump pulse to impul- 
sively align molecules By taking the HATI spec- 
tra at the time when the molecules are preferentially 
aligned or anti-aligned with respect to the polarization 
of the pump beam, the dependence of the DCS on the 
alignment of molecules can be obtained. Such measure- 
ments are beginning to emerge [55] ]. However, current 
interpretation of such experimental data, for molecules 
that are aligned [55| or not aligned (5(|, rely on intu- 



itive simple physical model 55 1, or on the extension of 
SFA2 to molecular targets 57, [H^]. As shown in the 
present paper, SFA2 is not expected to adequately de- 
scribe the backscattering of the returning electrons by 
the molecular ions. Using the QRS, more accurate de- 
scription of the HATI spectra is possible if the DCS from 
electron-molecular ions collisions are available. Electron- 
molecule as well as electron-molecular ion collisions have 
been studied over the past few decades, both in exper- 
iments and in theory. A few general purpose codes are 
in existence [H^, [H^, [62^ . By adopting these codes 
to obtain the DCS needed for the QRS model, accurate 
HATI spectra from aligned molecules can be calculated. 

Another potentially very important application is to 
use infrared lasers for dynamic imaging of a transient 
molecule [(33| ■ Using a pump pulse to initiate a chemical 
reaction, the time evolution of the transient molecule, 
including the position of its constituent atoms, can be 
probed by measuring the HATI spectra by a probe laser. 
Today few-cycle pulses of duration of a few femtoseconds 
are readily available, and as indicated in Section IV, no 
CEP stabilization is needed in order to extract the DCS. 
From the dependence of the DCS with respect to the time 
delay, one may be able to extract the time evolution of 
the structure of the transient molecule, thus achieving 
dynamic chemical imaging of molecules with temporal 
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resolution of a few femtoseconds. The feasibility of ex- 
tracting the structure of the target from the DCS has 
already been demonstrated for rare gas atoms using a 
simple genetic algorithm [32| . Extension of the method 
to molecular targets is likely to be straightforward. Thus 
while the idea of using laser induced electron diffraction 
for imaging the structure of molecules had been pro- 
posed since 1996 jfjij], due to the relatively low energies 
of the retur ning electrons, the standard electron diffrac- 
tion theory [65[ is not applicable. On the other hand, 
with the QRS proposed here, and with the implementa- 
tion of the state-of-the-art electron-molecule scattering 
codes, we believe that the fundamental theoretical tools 
needed for retrieving the structure of transient molecules 
based on laser-induced high-energy photoelectron spectra 
have been established. 



Before closing, perhaps it is fitting to mention that an 
expression similar to Eq. (1) has been used to explain 
high-energy electron emission in the forward directions 
in energetic ion-atom collisions, such as the collision be- 
tween Cu 5+ with H2 in Liao et al. [6^|. These electrons 
are best understood in the projectile frame where a beam 
of "free" electrons from H2 are incident on the Cu 5+ ion. 
When these electrons are backscattered by the ion, they 
emerge as high-energy electrons in the laboratory frame. 
In this model, the beam of "free" electrons from H2 is 
treated as a wave packet, represented by the Compton 
profile of the target, in the projectile frame, similar to 
the returning wave packet W(k r ) in the present QRS 
model. The transformation of momentum from the pro- 
jectile frame to the laboratory frame in ion-atom collision 
is equivalent to the — A r term in the present paper, see 
Eqs. (j48| and (j49|) . In fact the similarity is so close that 
we quote the abstract of here: "We present a method 
of deriving energy and angle-dependent electron-ion elas- 
tic scattering cross sections from doubly differential cross 
sections for electron emission in ion-atom collisions. By 
analyzing the laboratory frame binary encounter electron 
production cross sections in energetic ion-atom collisions, 
we derive projectile frame differential cross sections for 
electrons elastically scattered from highly charged projec- 
tile ions in the range between 60° and 180°. The elastic 
scattering cross sections are observed to deviate strongly 
from the Rutherford cross sections for electron scatter- 
ing from bare nuclei. They exhibit strong Ramsauer- 
Townsend electron diffraction in the angular distribu- 
tion of elastically scattered electrons, providing evidence 
for the strong role of screening played in the collision." 
Indeed, electron scattering experiments can be carried 
out without directly preparing a well-collimated electron 
beam! 
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